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ABSTRACT 

Recent progress in the measurement of relative distances to galaxies has 
been quite substantial, and catalogs of 3000 galaxies with distances are soon 
to become available. The peculiar velocity field (deviations from Hubble flow) 
derivable from these catalogs, when compared to the peculiar gravity fleld derived 
from all sky redshift surveys of galaxies such as the 1.2Jy IRAS survey, leads to a 
unique and extremely powerful test of the density parameter /3 = ^^'^ /hj, where 
hi is the possible linear bias of the IRAS selected galaxies relative to the mass 
fluctuations. We review the status of these large scale flow measurements and 
present a new methodology to describe the two flelds by means of an expansion 
in a set of orthogonalized functions describing a general potential flow to any 
chosen resolution. The parameters of the flow can be estimated by minimization 
of the describing the scatter of observed versus predicted linewidths from an 
inverse TuUy- Fisher relation. By this method one can intercompare the gravity 
and velocity flelds coefficient by coefficient, deriving a precise fit for the density 
parameter and an assessment of the degree of coherence between the fields. 
The present situation is transitory- different analyses of the same data are not 
yielding consistent results. Until this embarassment is untangled, estimates of 
13 should be taken with a large grain of salt. 

1. INTRODUCTION 

The study of large scale fiows permits an estimate of the density parameter 
on a characteristic scale of 10 — 30/?.~-^ Mpc, which is thirty times larger than 
the central regions of rich clusters, where mass estimates can be derived on the 
basis of hydrostatic equilibrium or virial equilibrium. As is well known, mass 
estimates from cluster analysis imply a cosmological density of order O 0.1 — 
0.2, assuming that galaxies in and out of clusters have similar mass to light ratios, 
or that the baryon fraction in the clusters is typical of the Universe as a whole 
(White eA al. 1993). The large scale fiow analysis would be expected to yield 
similar estimates of the density parameter if the dark matter fully congregates 
with the galaxies on cluster scales, but could well yield a different value if galaxy 
formation is more efficient in the denser regions or if the dark matter is too hot 
to cluster on any but the very largest scales. 

Dekel (1994) and Strauss and Willick (1995) have recently reviewed this 
subject in detail, and we shall confine our remarks here to a few comments on 
the discrepancies of the 13 derived from different analyses. A decade ago the 
quantitative analysis of large scale fiows consisted largely of the question of Vir- 



gocentric flow (Davis and Peebles 1983), but now the focus has broadened to 
analysis of the details of the flow on a scale four times as large, which has led 
to new insights about large scale structure in the Universe. The observational 
situation is improving rapidly, and will continue to improve in the next few 
years, particularly if high precision distance estimators, such as SN-Ia (Riess 
et al. 1994), or surface-brightness fluctuations (Tonry eA al. 1994) in early type 
galaxies, lead to useful catalogs. Giovanelli eA al. (1994) are preparing a full 
sky sample of some 2400 galaxies with TuUy-Fisher distances, while Willick eA 
al. (1995) have compiled a sample of 2900 galaxies over the full sky (Faber eA 
al. 1994). These large samples probe the velocity fleld to a redshift of approx- 
imately 8000 km/s, and the distance accuracy of the TuUy- Fisher indicator is 
estimated to have a precision of 15-20%. 

The cosmological interpretation of large scale flows is dependent on late 
time linear perturbation theory, which extensive testing has shown to be valid 
on scales > lQh~^ Mpc. This is convenient, because linear theory is easy to 
work with and is the fundamental prediction of alternative models of large scale 
structure. If one conflnes the analysis to simple moments of the velocity fleld, 
such as a bulk flow or the alignment of the flow on a given scale with the dipole of 
the cosmic microwave background, the cosmological test is limited to questions 
of the likelihood that a given model would produce sufficient power to generate 
such a flow. 

To sharpen the cosmological probe, it is instructive not simply to analyse 
the power in the velocity flow, but to compare the alignment of the velocity fleld 
with the inferred gravitational fleld derived from the observed galaxy distribu- 
tion. Toward this end, complete redshift surveys of whole-sky galaxy catalogs are 
an essential ingredient. The IRAS selected 1.2Jy survey (Fisher et al. 1995) of 
5300 galaxies covering 88% of the sky is the largest available catalog at present, 
but it should soon be superceded by at least two other catalogs (Santiago et 
al. 1994) (Saunders et al. 1995). The IRAS catalog is very suitable for this 
type of analysis because the 60 micron flux measured by IRAS is unaffected 
by galactic extinction and it is possible to select galaxies down to galactic lati- 
tude l&l = 5° with only modest confusion. However, the IRAS sample is quite 
dilute, particularly in cluster centers, and the inferred gravitational fleld has 
considerable noise, especially for structures with redshift cz > 6000 km/s. 

If one had full knowledge of the mass fluctuation fleld Sp(r) over all space, 
one could trivially write the gravitational fluctuation fleld g(r) as 



where p is the mean mass density of the Universe. The great simpliflcation of 
late time linear theory is the intuitive relation between the gravity fleld g and 
the velocity fleld vl, namely the equating of the velocity as the acceleration 
times the time, where the only possible time is the Hubble time. The correct 
expression (Peebles 1980) is 
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Assuming the galaxy distribution at least approximately traces the mass on 
large scale, with linear bias b between the galaxy fluctuations Sg and the mass 



fluctuations, (b = Sc/Sp), and replacing the integral over space with a sum over 
the galaxies in a catalog, with radial selection function </>(/•), we have 
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where n is the mean galaxy density in the sample. Note that the result is 
insensitive to the value of ifo, as the right hand side has units of velocity. Thus 
in this game it is traditional to quote distances in units of km/s. 

Because the effective weight per galaxy l/(/)(r) diverges at large distance, 
one must cutoff the calculation at some limiting redshift, which for the 1.2Jy 
IRAS sample is usually chosen to be 20,000 km/s. To limit the direct influence 
of one galaxy on another and to fllter out nonlinear effects, it is necessary to 
smooth the small scale interaction, typically with a smoothing scale of Vg > 500 
km/s. The method assumes that the associated mass per galaxy is constant, 
but in fact it is straightforward to piece together different redshift catalogs with 
different weights and selection functions. 

Two different schemes have been used to compute vl from redshift cata- 
logs. The original scheme of Yahil eA al. (1991) solved Equation (3) by iteration, 
adiabatically turning on the gravity fleld to advance the position of the galaxies 
from redshift space to real space. Recent reflnements on this method by com- 
putation on a grid (Strauss and Yahil 1994) will allow nonlinear corrections to 
be applied to the predicted velocity fleld. An alternative formulation that can 
be applied directly in redshift space has recently been suggested by Nusser and 
Davis (1994a). This method solves a modifled form of the Poisson equation using 
spherical harmonic decomposition and can readily deal with variable smoothing. 
To lowest order the two methods agree, but differ on how they treat the triple 
valued zones, regions around cluster centers in which three separate distances 
can have the same redshift. The Nusser and Davis algorithm intrinsically as- 
sumes a monotonic relation between distance and redshift, thereby eliminating 
multivalued zones. The older, iterative algorithm allows the existence of mul- 
tivalued zones but has no idea of how to position an object that flnds itself in 
such a zone; several ad hoc algorithms have been adopted over the years. There 
is no good solution for this problem which is one reason to fllter out nonlinear 
effects or to compare alternative algorithms. Fortunately, most of the local den- 
sity fleld appears to be single valued when smoothed on scales of 500 km/s or 
larger; only around cluster centers is the reconstruction of questionable validity. 
Nonlinear reconstruction algorithms are of marginal utility because they do not 
eliminate the problems of the multi-valued zones, yet the nonlinear effects are 
only important in regions which are nearly identical to the multivalued zones. 
Furthermore these are the regions where the assumption of linear bias is most 
questionable. 

2. S-S COMPARISON 

The POTENT algorithm (Dekel et al. 1990) is a scheme to convert the 
observed radial velocity fleld into a density fleld, and depends on the reasonable 
assumptions that the large scale velocity fleld is a potential flow in the linear 
limit. The assumption of potential flow allows one to construct a three dimen- 
sional velocity fleld from the observed radial flow, and the divergence of this fleld 
yields an estimate of the fluctuating mass fleld driving the flows. Comparison of 



this mass field with the observed galaxy distribution leads to an estimate of (3. 
Details are given by Dekel (1994); the method has been applied to the Mark-2 
catalog (Faber and Burstein 1988) by Bertschinger et al. (1990) and then com- 
pared to the galaxy density field derived from the 1.9Jy IRAS redshift catalog 

Strauss eA al. 1992), resulting in a very high density estimate, 13 = 1.3 it 0.3 

Dekel eA al. 1993). 

In order to define the velocity potential, it is necessary to heavily smooth 
the velocity data; Dekel eA al. (1993) use a guassian smoothing of cr = 1200 
km/s, which is sufficiently large to filter out the infall of the local group toward 
the Virgo supercluster! In its present implementation the POTENT method 
depends on a forward TuUy- Fisher approach, in which a calibrated TuUy- Fisher 
relation is used to infer the absolute magnitude of a galaxy given an observed line 
width. As is well known, this procedure suffers from homogeneous, and more 
particularly, inhomogeneous Malmquist bias (Schechter 1980). If the scatter 
in the distance estimator is accurately known, the Malmquist biases can be 
removed. If the incorrect scatter is used, of if no correction is applied, an 
overdensity of points in redshift space will appear to have the signature of a 
massive cluster, even if the cluster has no mass at all. Dekel eA al. (1993) 
undertook an extensive series of Monte-Carlo tests to calibrate these biases, 
and their result is consistent with a number of earlier comparisons of IRAS 
predictions to peculiar velocity data (e.g. Kaiser eA al. 1991, Straus 1988, Strauss 
and Davis 1988). 

3. v-v COMPARISONS 

Because this result is the only observational evidence for a high density 
Universe, it is certainly worth pursuing alternative methods that could provide 
confirming or refuting evidence. In the interval since Dekel eA al. (1993), the 
available database has improved considerably. The IRAS sample has doubled 
with the availability of the 1.2Jy survey (Fisher 1992), but more importantly, 
the Mark-2 sample of 496 galaxies is being superceded by the Mark-3 sample of 
2900 galaxies (Faber eA al. 1994, Willick etal 1995a). Several new analyses are 
underway and results will be reported in due course. Here we give a progress 
report. 

Willick eA al. (1995b) compare the observed distribution of the magni- 
tudes, redshifts, and linewidths for the Mark-3 data, and have cast the problem 
in terms of a likelihood of measuring a set of observed magnitudes, given the 
positions of the galaxies in redshift space and observed linewidths. The method 
accounts for the inhomogeneous space distribution by means of the observed low 
resolution density field measured from the IRAS sample. They have devised an 
algorithm that explicitly attempts to neutalize the effects of the triple valued 
zones in a statistical sense. Distortions from real to redshift space are provided 
by flow models derived from the IRAS gravity fleld, with 13 the only adjustable 
parameter. The procedure then asks which value of 13 is most likely to flt the 
observed distribution. The test favors a tentative value of /3 = 0.55 it 0.13. 
Details are given by Strauss and Willick (1995). While the method provides a 
powerful measure of which IRAS flow model is the best flt to the data, it does 
not in itself provide an independent picture of the flow fleld, which POTENT 
does provide, and which gives POTENT such a visual appeal. 

Two alternative schemes of intercomparing the velocity and gravity flelds 
have recently been presented by Nusser and Davis (1994a,b). The flrst is limited 
to a consideration of the dipole component of the radial peculiar velocity fleld 



and depends on a simple, but little used fact of potential theory. Consider a de- 
composition of the radial field on a given redshift shell into spherical harmonics. 
Let the shell have radius r and position a charge (mass) at some location R. 
For a potential field $ satisfying the Laplace equation, the l^^ multipole $; will 
scale as r' or /-"C+i)^ depending on whether R > r ov r < R, and the gradient 
of the potential (v or g) will scale as r^~^ or as 7'~('+2). "\Ye are used to the iron 
sphere theorem of Newton: for / = 0, a particle feels no force from the mass 
external to it. Note that for 1 = 1, d^i/dr is constant for r < R. That is, by 
transforming to the local group frame, the dipole amplitude of the velocity on all 
shells with r < i? is zero, and the influence of an external charge is measurable 
only for / > 2. In the local frame, the measured dipole on a given shell must be 
due entirely to the distribution of charge (mass) internal to that shell. 

Such a situation overcomes a serious objection to the use of Equation (3)- 
that the gravity fleld is nonlocal. Since the estimates derived from IRAS do not 
integrate the mass density over all space, they might be missing some large scale 
components. This problem is explicitly avoided by limiting the consideration to 
the behavior of the dipole amplitude. Nusser and Davis (l994a) extracted the 
dipole component of the velocity fleld from an early version of the Mark-3 PO- 
TENT map and compared it to the dipole fleld of the IRAS gravity maps, for 
various values of (3 (Figure 1). They flnd that the velocity and dipole vector 
directions are aligned within 11°, and furthermore that wiggles in the ampli- 
tude of the dipole versus redshift z are remarkably consistent. Note that the 
POTENT dipole does not go to zero at the origin because of the large smooth- 
ing. (3 = 0.6 lb 0.2 seems to be the preferred flt, in agreement with Willick et 
al. (1995). One very important point from Figure 1 is that 450 km/s of shear is 
clearly generated within a distance of 5000 km/s, a distance to which the IRAS 
density fleld should be reliable. Only the 250 km/s difference of this measured 
shear amplitude and the 620 km/s dipole amplitude of the CMBR dipole can 
be induced by larger scale bulk flows. This must be kept in mind in the inter- 
pretation of claims of the detection of a bulk flow on a much larger scale (e.g. 
Lauer & Postman 1994). 

4. INVERSE TULLY-FISHER RECONSTRUCTION OF A SMOOTHED FIELD 

To generalize this approach, one might consider directly comparing the 
velocity fleld derived from POTENT with the IRAS gravity fleld for higher 
multipoles. This approach is not advisable, because the distribution of points 
comprising the Mark-3 sample are not uniformly distributed on the sky or in 
redshift, so that the multipoles are not orthogonal. Furthermore POTENT is 
constructed on a cubic grid, which is not suitable for higher order multipole 
comparison. 

Consider an alternative approach that starts afresh with the measured 
magnitudes and linewidths. Nusser and Davis (1994b) describe an inverse TuUy- 
Fisher algorithm that based on a flow model which is a general, smoothed de- 
scription of the velocity fleld in redshift space. Assume that the TuUy Fisher 
relationship is linear, i.e. 

rii = ln{Avi) = 7]o + sMi ± e , (4) 

where 770 and s are the unknown zero point and slope of the regression, e is the 
error in the relation with (e) = 0, and 

Mi = nii — 5log(zi — Ui) — 15 . (5) 
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Fig. 1: The dipole amplitude of the Mark-3 sample compared to the IRAS pre- 
dicted dipole amplitude for various (3 values. 

Here mi is the apparent luminosity, Zi is the redshift, and Ui is the peculiar 
velocity of the i^^ object in the LG frame. Let Mi = Moi + Pi with Mio = 
mi — 5log(zi) — 15 and Pi = —5log(l — Ui/zi). Then, if the TF scatter of galaxies 
in the 77 direction is cr^, one can compute a statistic that minimizes the 77 
scatter as a function of the parameters that describe the peculiar velocity field 
P. Such an inverse fit to the TuUy-Fisher relationship does not suffer the usual 
Malmquist biases, but is biased in the presence of a "hot" velocity field or if 
there is a selection in the sample based on 77. The typical bias in the inferred 
velocity is Au^ias ~ ^m/^- Given that the rms scatter of distances at a given 
redshift (T„ is thought to be less than 300 km/s in most regions of redshift space, 
the bias of this technique is quite modest. 

The problem with inverse methods is that, in contrast to a forward method 
such as POTENT, one does not generate pictures, but can only fit parameters 
to a model. However, with large datasets such as Mark-3, which have nearly full 
sky coverage to substantial depths, this is not a serious limitation. Consider for 
example an expansion of P in terms of a set of functions that are orthogonal 
over a given dataset, 

i 

where the aj are the coefficients and F- is the value of the j^^ orthogonal 
function for object z, and 

i 




with Sj^ the Kronecker delta function. Orthogonal functions are a convenient 
choice, as the minimization is linear in all the coefficients aj as well as 770 
and s, and each term of the orthogonal expansion is decoupled from the others 
so that terms can be added one by one. 

A natural choice of functions for a general expansion of the radial peculiar 
velocity field is again motivated by potential theory. Because we are working 
in the local group frame, P(0) = 0. One can guarantee the correct asymptotic 
behavior near the origin by expanding the field P in terms of the derivatives of 
spherical Bessel functions times spherical harmonic functions Y;m(^5</')5 

nlm nlm 

See ND2 for details. One can choose kn by the desired boundary conditions at 
some czmax = 10000 km/s, such as requiring F(zmax) = 0, but the solutions 
are insensitive to the details of this choice since the Mark-3 sample sets little 
constraint on the flow at this redshift. Fisher et al. (1994) have used similar sets 
of wavefunctions to describe the density fleld of the IRAS 1.2 Jy catalog. The 
choice of radial wave function can be flexible, as the method is only meant to 
be a fltting formula, and the irregular distribution of galaxies in redshift space 
prevents the original basis functions from being orthonormal. A convenient 
choice that leads to lower with fewer degrees of freedom is to use radial 
wavefunctions j'(kny)/z where y = (ln(l + (z/2000)^)-^/^. This has the effect of 
stretching the oscillations of the radial wavefunctions at larger redshift, allowing 
more radial resolution in the foreground than in the background, so as to match 
the resolution gradient of the gravity fleld and the absolute accuracy gradient 
of the observed peculiar velocities. By working to order rimax = 5, Imax = 4, 
we have 125 independent modes. Five additional modes describing an external 
quadrupole fleld (F = constant X i^2m(^5 are also included for good measure. 

We transform the functions F- into an orthogonal basis set F- by means of an 
SVD inversion (Press et al. 1992). 

To conflrm that the method works, we constructed 2600-point mock Mark- 
3 catalogs derived from nbody simulations, including inhomogeneous sky selec- 
tion and scatter in the TF relation of cr^ = 0.03. Of the 130 modes, only 
approximately 80 are determined with signiflcance greater than one sigma, but 
these 80 coefflcients contain a full description of the velocity fleld to the chosen 
resolution. The scatter in 77 about the flt solution is spatially random in the 
simulations, indicating that flts to this resolution are adequate for the size and 
precision of the catalog. More detailed flts can be performed with larger set of 
TuUy-Fisher data, or with more precise distance indicators. 

To test whether the estimated mode amplitudes are correct, we have eval- 
uated each mode amplitude using the true peculiar velocity of each simulation 
point. Figure 2a shows the comparison of the true versus ITF inferred mode 
amplitude for the 130 coefflcients derived from the mock catalog. No bias in the 
mode amplitudes are detected, and the scatter of the amplitudes is completely 
consistent with the expected noise. Furthermore, most of the power of the large 
scale flow is contained within a few basis functions. For the true sky, we of 
course do not know a priori the true peculiar velocity of each point, and we 
can only substitute the linear-theory predicted peculiar velocity derived from 
the IRAS gravity maps. Figure 2b shows this test with the same mock Mark3 
catalog, now substituting an IRAS derived estimate of the peculiar velocity of 
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Fig. 2: (a) The ITF inferred amplitude of the orthogonal modes of a mock Mark-3 
catalog are plotted versus the true mode amplitude, (b) The solid symbols 
show the ITF inferred amplitudes versus the amplitude measured from the 
linear theory predictions for (3 = 1 (the correct answer in this simulation). 
The scatter is larger than Fig. 2a, but the slope remains unbiased. The 
open symbols are for a gravity field with (3 = 0.6; note that the diagonal 
line is now a poor fit. 

each point in the evaluation of the mode amplitudes. A flux limited mock IRAS 
catalog with the same selection function as the 1.2Jy survey was extracted from 
the simulation data and used to estimate the gravity fleld, so this test is a com- 
plete simulation of the real situation. Now a comparison of ITF inferred versus 
IRAS inferred mode amplitude shows more scatter, but still a complete absence 
of bias. The simulation is evaluated for (3 = 1, the true answer in this case. 
Evaluating the gravity fleld for other (3 values changes the slope of the scatter 
diagram, so estimation of (3 is completely straightforward. A gravity fleld evalu- 
ated using (3 = 0.6 is shown as the open symbols, and is clearly a poor flt to the 
diagonal hne. With a suitable catalog of distance measurements, this method 
should lead to a deflnitive measure of (3. 

The great advantage of this scheme is that it allows the comparison of 
the velocity and gravity fleld mode by mode, with identical resolution to each 
fleld. The Monte-Carlo simulations with realistic mock catalogs conflrm that 
the bias in the extracted velocity fleld is negligible. The method allows for the 
gradient in radial resolution of the two flelds, so it permits the extraction of the 
maximum useful information for each fleld. One can rank order the modes by 
their signflcance, or contribution to keeping only the most signiflcant modes, 
so as to provide a very compact description of the fleld. This is an example 
of a Karhunen-Lo'eve transformation (Therrien 1992). Recent applications of 
this technique to astrophysics include the CMBR fluctuations of the COBE data 
(Gorski 1994, Bunn 1994, Bond 1994), the power spectral analysis of fluctuations 
in redshift surveys (Vogeley 1994), and the spectral classiflcation of galaxies 
(Connolly et al. 1994). 

However, at the present time, this method is not yielding sensible results 
when applied to the Mark-3 catalog. There is strong coherence in the residuals 



between the ITF velocity field and the gravity field which is unphysical and 
suggestive of a calibration error. Work is in progress to iron out the problem 
and results will hopefully be available soon. 

5. WHY DON'T THE ALTERNATIVE METHODS YIELD A CONSISTENT (3 ? 

The POTENT-IRAS density comparison resulted in an estimate of /3 > 1, 
while the v-v comparisons discussed above are giving a best value half as large. 
What could account for such a discrepancy? 

It is important to note that, while the different methods might start with 
the same basic data, their complex processing is rather different. In the 8 — 8 
analyses, the weighting is dominated by infall to cluster centers and outfiow from 
voids. Larger scale components of the velocity field are less important than the 
local divergence. The direct velocity- velocity comparisons, on the other hand, 
are dominated by the signature of the dipole refiex of the motion of the local 
group. Figure 1 shows that 2/3 of the motion of the LG is refiected within the 
galaxy distribution to redshift 6000 km/s. Giovanelli eA al. (1994) are claiming 
the refiex dipole is larger than shown in Figure 1, which would increase the 
inferred 13 toward the POTENT value. The coherent residuals between the ITF 
inferred field and the IRAS predicted field for the Mark-3 sample are certainly 
cause for concern. Until we understand better the reason for this discrepancy, 
we must beware of possible systematic errors in 13 values so derived. 

The subject of large scale fiows is rapidly evolving. The large datasets from 
Giovanelli eA al. (1994) and Mark-3 are likely to be released soon, and other, 
more accurate catalogs of peculiar velocities (e.g. using surface-brightness fiuctu- 
ations) are under construction. The 1.2Jy IRAS catalog will be soon superceded 
by the IRAS catalog of 15,000 galaxies with fiux greater than 0.6Jy (Saunders 
eA al. 1994) and by the ORS/IRAS catalog (Santiago eA al. 1994), both of which 
will offer denser sampling and therefore improved statistical precision to the 
gravity field. The large scale fiow analyses should yield a 13 value accurate to 
10%, and the various estimates should certainly be consistent. Since progress 
continues at a rapid pace, it is perhaps premature to expect consistency on this 
important parameter at the present moment. 
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